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Abstract 

Diffusion of transported particles in the intracellular medium is described by means 
of a generalized diffusion equation containing forces due to the cytoskeleton net- 
work and to the protein motors. We find that the enhanced diffusion observed in 
experiments depends on the nature of the force exerted by the protein motors and 
on parameters characterizing the intracellular medium which is described in terms 
of a generalized Debye spectrum for the noise density of states. 
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1 Introduction 



The study of the spatio-temporal organization of vesicles, organelles and other 
particles inside the cellular medium plays a primary role in the behaviour of 
cells and constitutes a central problem in physical biology [US]- 

In order to carry out this study, it is important to consider the fact that 
the intracellular medium consists of a wide variety of polymers and particles 
whose presence is the cause of the observed viscolestic behaviour [3j. As a 
consequence, during their motion inside the cell the particles may experience 
elastic, electrochemical and entropic forces that introduce memory effects. 



E-mail: isholek.fc@gmail.com 



Preprint submitted to Physica A 



4 February 2009 



These effects must be taken into account in order to correctly describe the 
dynamics of the particles. 



Particles in the intracellular medium undergo different forms of diffusion which 
depends on the activity of the protein motors and the size of the particles 



[^ll5]|6H71l8|9|10llll|12|13|14|15|16j . Subdiffusion has been observed in the passive 

transport of particles in the cellular medium [H[5|9l[T0fTT] whereas an enhanced 
diffusion has been reported for particles driven by protein motors piip^ . 

A complete description of particle transport must take into account the ac- 
tivity of protein motors [3PfT2] . These entities are chemo-mechanical trans- 
ductors which, in the presence of sufficiently high cation concentration, trap 
particles and then use the free energy arising from the hydrolysis of ATP 
molecules to drive them along defined trajectories inside the cell pill7fl8j giv- 
ing rise to the enhanced diffusion observed in these systems [Sfl] . Experiments 
seem to justify the assumption that the activation as well as the duration of 
the activity of these protein motors constitutes a random process [3] . 

In this article, we use mesoscopic nonequilibrium thermodynamics P^2UII21] 
to formulate a mesoscopic model that accounts for the enhanced diffusion of 
particles in the intracellular medium. We propose a generalized diffusion equa- 
tion, that incorporates the memory effects due to the viscoelasticity of the 
intracellular medium and the activity of protein motors, to compute the mean 
square displacement (MSD) and to analyze the density of states of the intracel- 
lular medium [5ll9]|22] . Memory effects can be incorporated through effective 
mobility and diffusion coefficients that can be memory kernels, [5|22|23ll24j 




or time- dependent coefficients [6|f7|l8|l9|l22] . In our approach, the generalized 
diffusion equation incorporates memory effects through time dependent coef- 
ficients. We choose this way because in the linear force case it is possible to 
show the equivalence of this description with that given by means of gener- 
alized Langevin equations with memory kernels [6|7]. This relation between 
both descriptions is also used in order to calculate the noise density of states 
arising in the case of enhanced diffusion. However, in this case the stochas- 
tic forces are of non-thermal nature due to the activity of protein motors. It 
is convenient to mention here that, in the general case, there is not a clear 
criterion in order to choose between time-dependent coefficients and memory 
kernels, see, for example, Ref. [25] . 

The article is organized as follows. In Sec. 2, we derive a generalized diffusion 
equation (GDE) describing the dynamics of the system. In section 3, we an- 
alyze the case of a harmonic force model and compare the results obtained 
with experiments. Section 4 is devoted to derive the explicit form of the noise 
density of states (NDS). In Sec. 5 we summarize and discuss our main results. 
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2 Mesoscopic thermodynamics approach to anomalous diffusion 



We consider a Brownian particle moving in an intracellular medium having 
viscoelastic nature [5]. In addition, we will assume that under the appropriate 
conditions the particle can also be driven by protein motors. As we have men- 
tioned previously, these factors introduce memory effects that will be taken 
into account through the effective time dependent mobility and diffusion co- 
efficients. 

These effects can appropriately be described in terms of a generalized diffusion 
equation for the single particle distribution function f[f,t), which gives the 
probability of finding the particle at position r at time t. This distribution 
function satisfies the normalization condition / /(r, t)df= 1, and its evolution 
in time is governed by the continuity equation 

|:/(f,t) = -V-[/(f,t)K(f,t)], (1) 



which expresses probability conservation and where V = d/df. The quantity 
fVr represents an unknown probability current and the explicit expression 
of the streaming velocity Vr can be calculated with mesoscopic nonequilib- 
rium thermodynamics (MNET) using the generalized Gibbs entropy postulate 
I51IT91I20] 



Sit) = -kB I f In ^j^df+ So, (2) 
fo{r) 



where ks is Boltzmann's constant and fo{r) and 5*0 are the distribution func- 
tion and the entropy of the reference local equilibrium state. The local equi- 
librium distribution function is given by fo{f) = eit^o-<t>)im/kBT) ^ which m is 
mass of the particle, T is the temperature of the system, 0(r) is an external 
potential per unit mass and fiQ is the chemical potential per unit mass at local 
equilibrium . 

Taking the time derivative of Eq. ([2]) and using Eq. ([T]) one obtains the ex- 
pression dS(t)/dt = —ks /ln[///o] V ■ {fVr)dr. Now, integrating by parts 
assuming that the current vanishes at the boundary, we find the entropy pro- 
duction a{t) 

ait) = -jJ f{f,t)Vr{f,t) ■ V^{f,t)df> 0, (3) 

where /i(r, t) = (kBT/m) In ///o + /io is the nonequilibrium chemical potential 
per unit mass. According to nonequilibrium thermodynamics, we will assume 
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that the force V/x(r, t) and the current f{r,t)Vr{r,t) are coupled hnearly: 
fVr oc V/U, [5]|9|20II21] . Using the expression for /o(r) one obtains 



f{r,t)Vr{r,t) = (3-\t)f{f,t) (^F{r) 



m 



Vln/(f,t) 



(4) 



where F{f) = — V0(r) is the total force applied on the particle. As men- 
tioned previously, the time dependent Onsager coefficient P~^it) accounts for 
the presence of memory effects. In general, this quantity is related with the 
correlation function of the position vector of the particle. We have chosen this 
coupling between the current and force because in the linear force case our 
results are equivalent to those obtained from a GLE in Refs. [7,8]. The case 
when the coupling is given through a memory kernel has been discussed in 
Ref. [6]. In the absence of memory effects, P~^(t) reduces to the Stokes mobil- 
ity /3q^ = (Qiraris/m)"^ with a and m the radius and mass of the particle and 
ris the viscosity of the medium. p!26] Substituting Eq. (jlj) into ([T]), we finally 
obtain the generalized diffusion equation 

|/(f,t) = ^p-\t)V'fif,t)-P~\t)V-[fir,t)Fir)] . (5) 
at m L J 



A similar equation has been proposed to describe finite-size and confinement 
effects on the dynamics of passively diffusing particles moving through the 
intracellular medium. [5|9] 



3 Enhanced diffusion in the intracellular medium 



A physical model describing the dynamics of a particle driven by protein 
motors must take into account two forces. The first one is a trapping force 
of elastic nature which in first approximation can be modelled by the linear 
force: Feiif) = —(jJ^r, with uq a characteristic frequency. The second one is 
the driving force exerted by molecular motors on the particles that can also 
be modelled by Fpm{r) = a;^r, with a characteristic frequency ui. 

— * — * 

Substitution of the total force F{r) = Fei{r) + Fpm{r) into Eq. ([5]) yields 



d_ 
di 



f{r,t) = p-\t) 



{ul - ■ [ff{f, t)] + ^ VV(r, t) 

m 



(6) 



where we have assumed that the effective friction coefficient is given by (3{t) = 
l3ol3{t). P{t) is a dimensionless function accounting for memory effects. Eq. ([6]) 
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clearly manifests the competition between trapping and driving forces that 
will regulate the behavior of the MSD of the particle. By introducing the 
dimensionless variables x = a~^r and i = r^^t, with td a characteristic time, 
and scaling time with z{t) = (3~^{t')dt' , we obtain 



e-2K--!)^-{t) 



(7) 



where we have introduced the constant uo^j. 



kBT/ma? and {x'^ii)) = J x'^fdx. 
Eq- gives the behavior of the MSD in the case when uJq — uf > 0, that is, 
when the friction is larger than the force exerted by the protein motors. The 
case when Uq — uf < will be described later. 



The time dependence of z{t) is related to the behavior of the time correla- 
tion function = {x ■ Xo)(t) in the absence of molecular motors, in the 
subdiffusion regime ui = 0. The evolution equation for xii) is given by Pl7f5] 

ixW = -^o^r^WxW, (8) 



which can also be calculated by using the dimensionless form of Eq. ([6]). Now 
by taking into account the expression of z{t) in terms of after solving ([8]) 
one arrives at the relation 

Zsubii) = -^\nR{i), (9) 



where we have defined the normalized correlation R{t) = x(t)/x(0). From this 
relation it follows that the scaling of time depends on the relaxation dynamics 
of the system. 

At short times, one may obtain the behavior of z{t) by expanding the logarithm 
in its argument around one: In |i?-*/^o'^o | ~ R-f^o/roi^o - 1 + 0{R'^). Then 
taking into account that R{t) must be an even function of time [U], a first 
order expansion leads to: R~^ii) ~ 1 + + O(t^). Thus, at short times we 
obtain the expression 

^.n6(t)~i?2p''°/"-'^°, (10) 



where B2 oc -B2 is a parameter characterizing the type of relaxation [9j. Now, 
by expanding the exponential in Eq.© up to first order in its argument and 
substituting (fTU]) into the result, we find 

(x2(t)).„,^652rD^P*/-^o. (11) 
Po 
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The exponent characterizing the subdiffusion process of the particle is a func- 
tion of parameters of the bath and the forces acting on the particle pTIlQ] . 
The coefficient B2 depends in general on the concentration of polymers in the 
medium and the size of the particles, thus characterizing the magnitude of the 
MSD, [9]. For t = Td, Eq. ([II]) yields B2 = {l3o/rpuj^o){^\'^))jub. }Jsmg the 
definition of the running diffusion coefficient Dsubii) = d{x^){i)/di, [22j, Eq. 
([7]) leads to the relation 

/3-Ut)^ri+^'3»/--o. (12) 



This expression can be used to obtain the shear modulus characterizing the 
viscoelastic properties of the material [5,..9]. 

In Fig. 1, we compare our results with experiments (symbols) reported in Ref. 
[T|, where the motion of a microsphere through a living eukaryotic cell was 
observed with video-based methods. The value of uq can be determined by 
considering the case when the particle performs subdiffusion or is attached 
to the cytoskeleton {uji = 0). From the experimental results represented by 
open circles in Fig. 1 it follows that the short time behavior of the MSD 
satisfies a power law with exponent ~ 3/4, [^10|27] . Using the parameters B2 
and td to fit the curve (dot-dashed line) with Eqs. ([7]) and ffTOj) . one obtains 
uJq ~ (8/3)/3or^\ implying that luq is proportional to the geometric mean 
of the fast ~ Po and slow ~ r^^ modes (see the caption of Fig. 1). Since 



Uq describes subdiffusion, it must satisfy the condition uq > y2(3o/TD. It is 
worth noticing that the exponent depends on the friction coefficient /5o, the 
characteristic frequency uq of the elastic force exerted by the medium on the 
particle and a scaling time td which is related with the time at which the 
MSD saturates. In a previous study this power law behavior was obtained by 
simply inferring the exponent from direct comparison with the experimental 
results Here, we show that it is connected with the interactions between 
the particles and the medium. 

The observed enhancement of diffusion can be described by taking into account 
the activity of the motors {uf 7^ 0) in Eq. IQ. Defining u'^ = — > 0, 
we use Eq. (|6]) in order to obtain a linear evolution equation for {x{t)) = 



J xf{x,t)dx whose solution is {x{t)) = exp oj'^tdISq ^ J Psuti'^')^^' 



Assuming 



that u'^tdPo ^ < 1, an expansion of the exponential up to first order yields 
(f(f)) ~ 1 + Cu^TDl3^'{x^ii))sub, (13) 



where we have used the relation Zsub(t) = J Psubi^)^^- Considering that, in gen- 
eral, protein motors are randomly distributed in the cell, we will assume that 
during the motion of the particle the parameter uj"^ is a fiuctuating quantity 
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with zero mean and {uj'^{t)uj'^{t')) = (exp —{t — t') , with the correla- 
tion time of the noise and ( = uf — its magnitude. Now, by evaluating 
{x{t)) {x{t')) and taking the average over the realizations of the random force, 
represented by {x{t))'^, one finally obtains 



(14) 



Fig. 1 shows a comparison between experiments (filled symbols) and theory 
(solid and dashed lines) of enhanced diffusion due to the presence of protein 
motors. For 2Po/ (uqTd) = 3/4, from Eq. IHM it follows that the exponent 
characterizing the enhanced diffusion is a^d ~ 3/2. The results ffTTj) and f|T^ 
constitute an alternative way to calculate the Hurst exponent of the MSD typ- 
ically arising in systems in which self-avoiding random walks (SAWR) appear 
frequently, such as crystallization or cluster-cluster aggregation in a viscoelas- 
tic medium [5|28f29ll30j . This fact allows one to interpret the present results 
within the scope of SAWR formalism [28j in which the exponent obtained is 
3/2 = 2/d with d the dimension of the space. The value of d corresponding 
to our problem in the enhanced diffusion case is therefore d = 4/3. Thus, the 
molecular motors act as elastic scatters inducing a SARW in a space with 
slightly reduced (motion-oriented) degrees of freedom. 



4 The noise density of states of the intracellular medium 



A global characterization of the intracellular medium can be performed through 
the noise density if states (NDS). This quantity contains in average form the 
statistical properties of the intracellular medium, assimilated to an effective 
medium in which both thermal and non-thermal (motor-induced) fluctuations 
control the dynamics of the transported particles. 

The NDS can be analyzed by means of a generalized Langevin equation p^llTS] 
which is equivalent to the following generahzed diffusion equation piS] 



-/(f,t) = Z},^;(t)VV(r,t), (15) 



where Deff{t) = d{x{t))'^ /dt is an effective diffusion coefficient incorporating 
the effects of the motors in average form. Eq. f|T5l) is equivalent to the general- 
ized Langevin equation for the velocity v{t) of the particle in the overdamped 
case 08] 

t 

J U{i - t')vit')dt' = F^(t), (16) 
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where we have considered the one-dimensional case for simphcity and is a 
stochastic force satisfying the conditions {F^) = 0, {F^ x{0)) = and 

{F^{i)F^{i')) = {v\,U{i-i% (17) 



with (f^)eg = lii[n.t^oo{v^)(t)- By Laplace transforming Eqs. f|T5|) and f|T6l) 
we can obtain the relation between the diffusion coefficient and the memory 
function in the following form 



Deffit) 



m 



£-1 



.sfl(s). 



(18) 



where n(s) is the Laplace transform of n(t) and C ^[fl{s)] the inverse trans- 
form. Using now the definition of Deff{t) in terms of Eqs. f|T^ one obtains 



This equation can be used to calculate the NDS of states p{uj) of the intracel- 
lular medium by decomposing the stochastic force F^{t) of the bath in a set 
of harmonic oscillators p^llTS] . The relation between the memory function and 
the NDS is found through the fluctuation dissipation theorem. One obtains 

U{i) = [ p{u)cos{ut)duj, (20) 



where uj is now a dimensionless frequency. 

Taking the Fourier cosine transform of Eq. fl2UI) with 0191] . the resulting ex- 
pression for the NDS of the intracellular medium is 

^|-l+4/3o/^g^D. (21) 



p{uj) ~ — r 

TT 



1 + 



4/^0 



sm 



From this equation it follows that the exponent controlling the behavior of 
the NDS depends on the frequency characterizing the elastic forces {ujq) in the 
intracellular medium. It is important to stress that the power law dependence 
on the frequency of the NDS given by fl2Tl) and valid up to the correspond- 
ing Debye cutoff frequency, is of the form of the generalized Debye spectrum 
proposed in Ref. [I^. 

The results we have obtained enable us to conclude that the activity of protein 
motors in the intracellular medium modifies the dependence on frequency of 
the NDS. For example, in accordance with experiments [SP] the MSD of the 
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particles grows by following the power law ~ implying that the 

exponent of the generalized Debye spectrum is 1/2, which contrast with the 
spectrum for subdiffusion which is —1/4. 



5 Conclusion 

In this paper, we have introduced a general formalism offering a multiscale 
description of transport processes in an intracellular medium. By taking into 
account the forces exerted the cystoskeleton and the protein motors in a gen- 
eralized diffusion equation, our model explains both enhanced diffusion and 
subdiffusion observed in experiments in living cells. The exponents associated 
to these behaviors are given in terms of the parameters characterizing the sys- 
tem. These results could be interpreted within a more general context in which 
the particle is transported through the elastic medium undergoing a random 
path through elastic scatters that induce a self-avoiding random walk |28j . 
This more general scenario permits to relate the diffusion in the intracellular 
medium studied with other processes in which self-avoiding random walks are 
an important ingredient. 

Finally, we have proposed the density of states as a useful quantity to describe 
the state of the cell. In the presence of protein motors, this quantity behaves 
as the Debye spectrum with an exponent larger than one. The model proposed 
may describe and quantify the observed different forms of diffusion of particles 
in a intracellular medium. 
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Fig. 1 Mean square displacements of microspheres engulfed in living eu- 
karyotic cells. Symbols represent data obtained with the fitting 
formula (x^) = obtained from experimental measurements 
in Ref. [1]. Circles correspond to the case of absence of motors 
with b = 0.026 and a = 0.75. Triangles and squares correspond 
to the case of presence of motors with b = 0.22, a = 1.5 and 
b = 0.009, a = 1.42, respectively. The dot-dashed line (subdifFu- 
sion) is the theoretical result given by Eq. ffTTj) using ffTOj) for the 
following values of the parameters: Dq ~ ksT /m(3Q ~ 10~^yum^s~^, 
/3o ^ 10^s-\ m ~ 10-^2^, a ~ 1/im, ujI ~ 3.4 ■ Ws'^, B2 ~ 4.36 
and To — 2.67s. The dashed and solid lines show enhanced diffusion 
obtained with Eq. ([H]) for ui ~ 4.2 ■ lO^s"^ and ui ~ 3.4 ■ lO^s"^ 
respectively. 
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Figure 1. 
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